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Tracking the momentum flux of a CME and quantifying its 
influence on geomagnetically induced currents at Earth 
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Abstract. We investigate a CME propagating towards Earth on 29 March 2011. This 
event is specifically chosen for its predominately northward directed magnetic field, so 
that the influence from the momentum flux onto Earth can be isolated. We focus our 
study on understanding how a small Earth-directed segment propagates. Mass images 
are created from the white-light cameras onboard STEREO which are also converted into 
mass height-time maps (mass J- maps). The mass tracks on these J-maps correspond to 
the sheath region between the CME and its associated shock front as detected by in situ 
measurements at LI. A time-series of mass measurements from the STEREO COR-2A 
instrument are made along the Earth propagation direction. Qualitatively, this mass time- 
series shows a remarkable resemblance to the LI in situ density series. The in situ mea- 
surements are used as inputs into a 3D magnetospheric space weather simulation from 
CCMC. These simulations display a sudden compression of the magnetosphere from the 
large momentum flux at the leading edge of the CME and predictions are made for the 
time-derivative of the magnetic field (dB/dt) on the ground. The predicted dB/dt were 
then compared with observations from specific equatorially-located ground stations and 
show notable similarity. This study of the momentum of a CME from the Sun down to 
its influence on magnetic ground stations on Earth is presented as preliminary proof of 
concept, such that future attempts may try to use remote sensing to create density and 
velocity time-series as inputs to magnetospheric simulations. 



1. Introduction 

The influence of a southward orientated magnetic field of 
the solar wind impinging onto the Earths magnetosphere is 
known to be the main driver for coupling the solar wind to 
the terrestrial system [e.g. Gonzalez and Mozer, 1974; Bar- 
gatze et al., 1985; Tsurutani et al., 1992a]. In this case, the 
oppositely oriented magnetic fields between the solar wind 
and the magnetosphere are more easily able to reconnect and 
thereby transfer energy, mass and momentum to the Earth. 
For this reason significant emphasis on understanding and 
predicting the magnetic field orientation has been pursued 
to improve our space weather predictive capabilities. 

The role of the solar wind density, N sw , in space weather 
studies is more complex. Statistical studies aimed at cou- 
pling the solar wind density to geomagnetic indices such 
as Dst show a weak correlation [O'Brien and McPher- 
ron, 2000]. Others suggest that N 3W should be significant 
in mediating the energy transferred to the magnetosphere 
[Borovsky et al., 1998; Thomsen et al., 1998]. The den- 
sity has been demonstrated to affect the magnetospheric 
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response by way of the solar wind dynamic pressure [mo- 
mentum flux, e.g. Xie et al, 2008]. A sharp increase in so- 
lar wind dynamic pressure can rapidly compress the Earth's 
magnetopause causing a large time rate of change in the 
local ground-based magnetic field intensity. While this ef- 
fect is not often considered as a measure of 'geo- effective' 
events, the space weather consequences can be severe. An 
example of this can be found in the 4th August 1972 event, 
where the peak Dst index was moderate compared to result- 
ing outage of the AT&T telecommunications cable and the 
geosynchronous satellite solar cell damage [Anderson et al., 
1974; Lanzerotti, 1992; Tsurutani et al., 1992b] . 

When investigating the drivers of space weather at Earth 
and in particular geomagnetically induced currents (GICs), 
it is important to distinguish between the different solar 
wind structures that initiate the chain of events that may 
lead to significant socio-economic losses. Different interplan- 
etary structures have been previously reported to produce 
significantly varying responses in the Earth magnetosphere 
[e.g. Huttunen et al., 2008; Miyoshi and Kataoka, 2005; 
Borovsky and Denton, 2006; Denton et al., 2006; Lavraud 
and Borovsky, 2008], it is therefore reasonable to expect dif- 
ferent responses for GICs. The two main drivers of GICs are 
from intense magnetic storms associated with interplanetary 
coronal mass ejections (CMEs) and co-rotating interaction 
regions (CIRs). Kataoka and Pulkkinen [2008] and Borovsky 
and Denton [2006] both found that the biggest problems for 
ground-based conducting systems were from CME drivers, 
whereas CIR-driven storms had a relatively minor effect. 

Fast moving interplanetary CMEs can be considered to 
be made of two primary constituents: the sheath region be- 
tween the shock front and the leading edge of the CME; and 
the ejecta itself. The ejecta is often magnetically dominated 
(i.e. plasma /3 <1) and displays the magnetic properties of 
a flux rope (FR) when seen in situ [Burlaga, 1988]. The 
sheath region however contains denser compressed plasma 
with a higher value of plasma /3. The direction of the mag- 
netic fields within this region is variable, but is often found 
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to vary within a 2D plane [Nakagawa et al, 1989; Jones 
et al, 2002; Kataoka et al, 2005; Savant et al, 2011]. The 
momentum flux, Pd yn = N SW V 2 within the sheath is typi- 
cally high. Huttunen et al. [2008] concluded that as a geo- 
effective CME propagates over the Earth (usually for a pe- 
riod of ~ 24 hours), the most intense GIC activity recorded 
on the ground is most likely to occur at the beginning, dur- 
ing the passage of the sheath region, between the shock and 
the leading edge of the magnetic obstacle. 

A high solar wind density can cause a change in the 
compression ratio of the bow shock, as is frequent within 
a geomagnetic storm [Lopez et al, 2004]. For southward 
magnetic field this serves to increase the sensitivity of the 
predicted effects to N sw [see also Lavraud and Borovsky, 
2008]. Recent studies have considered the effect of precon- 
ditioning the magnetosphere to increase the space weather 
effects. For geomagnetic storms driven by both CMEs and 
CIRs, Lavraud et al. [2006] showed that the Dst were under- 
predicted by a model for intervals that were preceded by an 
extended period of northward IMF. This is consistent with 
the hypothesis that a preconditioned high density plasma 
sheet will cause a larger than expected geomagnetic storm, 
as also expected from simulations [Lavraud and Jordanova, 
2007]. The density of the Earth's plasma sheet is regulated 
by the solar wind density and lags behind by a few hours 
[Borovsky et al., 1998]. For CME-driven storms the plasma 
sheet is generally more dense and persists for longer than 
those driven by CIRs [Denton et al., 2006]. A magneto- 
hydrodynamic (MHD) simulation studied by Siscoe et al. 
[2002] found that the polar cap potential saturation follows 
a power law with Pa yn - This then led Xie et al. [2008] to 
develop a model where the magnetosphere can be precondi- 
tioned by Pdyn. 

The Dst index is often used as a measure of the severity 
of a geomagnetic storm. This index is an average estimate 
of the global response as determined by four equatorially 
located ground stations. The Dst measures the horizon- 
tal component of Earth's ground magnetic field and is pre- 
dominately monitoring the strength of the ring current in 
the magnetosphere [Fukushima and Kamide, 1973; Liemohn 
et al, 2001]. Since a denser and cooler plasma sheet pro- 
duces a stronger ring current, CME-driven storms often pro- 
duce more intense Dst values. Also, low latitude (< 56°) 
aurora, as measured in Japan were found to be predomi- 
nately CME-associated [Shiokawa et al, 2005]. During a 
CME-triggered geomagnetic storm, the Dst profile often be- 
gins with a short duration rise which is triggered by the 
arrival of the CME-shock front (sudden storm commence- 
ment, SSC). The duration of the rise phase between the 
SSC and storm onset (SO) often corresponds to the sheath 
region upstream of the magnetic ejecta of the CME. The 
main phase of the geomagnetic storm, which corresponds 
to the magnetic ejecta traversing over the magnetosphere, 
is where a large decrease in the Dst value occurs. However, 
the phases of a geomagnetic storm can be time shifted under 
conditions where a strong southward Bz persists within the 
sheath region. The amplitude of the Dst decrease during 
the main phase is closely related to the southward magnetic 
field in the IMF and its subsequent magnetic reconnection 
with the magnetosphere. 

The momentum flux of a CME, and in particular the 
sheath region, is higher than that of the ambient solar 
wind. This introduces a changing pressure onto the mag- 
netosphere, and in the case of a CME-driven shock front, 
the change is abrupt. The subsequent compression of the 
magnetosphere is expected to generate currents within the 
terrestrial system which are then detected on the ground as 
magnetic field fluctuations. As displayed by a Dst time pro- 
file during a geomagnetic storm, the geomagnetic field un- 
dergoes fast changes over a short time period. The largest 
time derivative of the ground magnetic field, dB/dt, occurs 



during the initial phase (between SSC and SO) or during 
substorms. As indicated by Faraday's law of induction, 
dB/dt is key in estimating the geomagnetically induced cur- 
rents in technological conductor systems like power grids 
[Pulkkinen et al, 2006; Committee On The Societal and 
Economic Impacts Of Severe Space Weather Events, 2008]. 
Therefore dB/dt is often considered as a reasonable proxy 
for localised GIC activity [Viljanen et al, 2001]. In order to 
calculate the geoelectric field, additional information on the 
ground conductivity and specifics on the technological sys- 
tems involved are required [Pirjola, 2002; Pulkkinen et al, 
2007]. In this paper we will focus on tracking the effects of 
CME momentum flux to the proxy of GICs, namely dB/dt. 

Significant geomagnetically induced currents have been 
observed to affect technological systems in mid- and low- 
latitude region [e.g. Ngwira et al, 2008; Liu et al, 2009; 
W atari et al, 2009]; and therefore it is becoming increas- 
ingly clear that these geomagnetically storm-driven GICs 
are not just a high latitude phenomenon [Pulkkinen et al, 
2010a]. The extent to which GICs can be expressed by 
dB/dt is also not as straightforward as first expected due 
to the complications of estimating the ground conductivity 
and how this varies at different depths. Pulkkinen et al. 
[2010a] showed that in Japan where the tectonic plates are 
geologically active, the local subduction zone is able to af- 
fect GICs in an unusual way by following the amplitude of 
the local geomagnetic field rather than the time derivative; 
this highlights the complexity of predicting GICs. 

In this paper we focus on better characterising the N sw 
from remote observations as a solar wind driver for dB/dt 
of the geomagnetic field. We take steps towards creating a 
time series of the solar wind density, and show that the mo- 
mentum flux of a CME is an important parameter when pre- 
dicting space weather incidents. We show this by comparing 
the ground station data with predicted estimates from the 
space weather modeling framework (SWMF) [Toth et al, 
2005]. 

2. Remote-sensing observations 

The STEREO mission, launched in 2006 [Kaiser et al, 
2008], consists of two spacecraft that follow a trajectory sim- 
ilar to that of the Earth. As they separate from each other 
at a rate of ~ 45° per year, one spacecraft travels ahead 
of the Earth (ST- A) while the other lags behind (ST-B). 
Each spacecraft carries the Sun Earth Connection Coronal 
and Heliospheric Investigation [Howard et al, 2008, SEC- 
CHI,] imaging package, which contains an Extreme Ultravi- 
olet Imager (EUVI), two coronagraphs (COR-1 and COR- 
2), and the Heliospheric Imager (HI). The HI instrument 
on each STEREO spacecraft is made up of two wide-field 
visible-light imagers, HI-1 and HI-2 [Eyles et al, 2009]. The 
fields of view of HI-1 and HI-2 are of 20° and 70° angular 
extent, respectively, and under ordinary operation are nom- 
inally centred at 13.7° and 53.4° elongation in the ecliptic 
plane. Thus the ecliptic plane corresponds to a horizontal 
line that runs through the centre of the fields of view. Figure 
1 displays the location of the STEREO spacecraft in rela- 
tion to the Sun and Earth on 27th March 2011. The shaded 
regions indicate the field of view for the HI-1 cameras on 
both spacecraft. The direction of propagation for the CME 
nose is estimated from Jmap techniques [e.g. Sheeley et al, 
1999; Rouillard et al, 2008; Davies et al, 2009; Savani et al, 
2009, 2012a] and detailed below in § 2.3. 

2.1. CME mass calculation 

In order to convert the white light images that con- 
tain photometric information per pixel into mass per pixel 
we exploit the Thomson scattering properties of electrons 
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from the corona and inner heliosphere. The total light ob- 
served by the coronagraphic and heliospheric imagers on- 
board STEREO are from photospheric photons scattered 
by all the electrons along the entire line of sight (LOS). 
The difference in contributions of each electron along the 
LOS depends on the distance and the Thomson scattering 
mechanism [Billings, 1966; Vourlidas and Howard, 2006], 
therefore an estimate of the total number can be generated 
from white-light images [Hayes et al, 2001]. In order to es- 
timate the mass of a CME, the brightness contribution of 
the transient must be isolated from the background coronal 
signal. Previous studies have achieved this by subtracting 
a suitably-chosen pre-event image from an image containing 
the CME [Stewart et al, 1974; Howard et al, 1985; Poland 
et al, 1981]. A comprehensive explanation of the proce- 
dures required to make estimates of a CME mass from the 
raw telemetry data received on the ground is given by Vourl- 
idas et al. [2010]. However for completeness, the main key 
points in the procedure are given below: 

1. First, the relevant time stamps for images contain- 
ing the CME are chosen as well as a pre-event image. The 
pre-event image is ideally the last possible image prior to 
the CME entering the field of view in order to make the 
appropriate corrections for background. It is important to 
minimise the effects from evolutionary changes and solar ro- 
tation, and not to include another CME or other transient 
effects in the pre-event. 

2. The pre-event image is subtracted from the sequence of 
CME-event images. These calibrated images now display in- 
formation on the excess (depletion) of light in units of mean 
solar brightness, MSB (i.e a base-difference sequence). 

3. The excess (depletion) of MSB within each pixel in 
the image is converted into a number of electrons by using 
the Thomson scattering equations and by assuming all the 
electrons are located on a single plane determined by the 
fixed-phi J-map technique. 

4. The mass per pixel is calculated from each image by 
assuming a solar wind distribution of 90% H and 10% He. 
This corresponds to a mass of 1.97 x 10~ 24 <? per electron 
[Hildner et al, 1975]. 

5. The mass of the CME is then estimated by summing 
up the values of all the pixels containing the CME. 

Figure 2 is a combined and cropped image from the ST-A 
and includes a frame from COR-2, HI-1 and HI-2. Practi- 
cally speaking, the pre-event image is usually subtracted 
from the frames which occur later in time (i.e. when the 
CME has entered the field of view). However, this frame 
may also be subtracted from images that occurred earlier. 
In figure 2, the image from the HI-2A camera is one that is 
actually earlier in time than the pre-event image. Therefore 
we are demonstrating the feasibility of using the pre-event 
image as a static background for frames earlier and later in 
time. 

The pre-event image subtraction process is most suitable 
for relatively short time-scales due to the steady state as- 
sumption. This assumption therefore begins to breakdown 
for the long duration of the CME within the large field of 
view. Also, HI-2 is sensitive to the background star field. 
This means that the ability to detect a CME motion is re- 
duced if the star field is not appropriately removed during 
the initial image processing prior to making mass calcula- 
tions. This paper is focusing on the potential to track the 
momentum flux and to monitor the relevant effects at Earth. 
For this reason we choose to use the simplest image process- 
ing at this stage in order to emphasise the minimum capa- 
bility for space weather forecasting. Further processing of 
the images should be able to improve the tracking of the 
CME momentum [Howard et al, 2012; Howard and DeFor- 
est, 2012] and may provide a more sensitive time series of 
mass flux at the LI point. 

The majority of the current work on mass estimates use 
data from coronagraphs which have a small field of view and 



assume the Thomson sphere is a flat plane. For the large 
field of view for both the HI cameras this is inappropriate. 
When estimating the mass, the propagation direction, and 
therefore the angle away from the plane of sky (PoS), is 
important in calculating the amount of Thomson scattering 
from electrons; i.e. the pixel on the inner edge of the camera 
has a different angle away from the PoS to the outer edge. 
This modification to the calculations is made in our work. 
In this paper, we treat each pixel individually as a different 
angle away from the plane of sky as measured along the line 
of sight. 

2.2. Semi-automated CME mass 

The calculation of the total CME mass is carried out by 
summing up the mass values of each pixel within the ob- 
served CME. This region can be defined in a variety of ways: 

1. The sector method [Vourlidas et al, 2010]. The ob- 
server manually defines a set of four boundaries that can 
be used for the entire sequence of images. These boundaries 
are defined between two Position angles (P.A.) and the inner 
and outer radial boundary. 

2. Region of interest (ROI) method [Vourlidas et al, 
2000; Subramanian and Vourlidas, 2007]. The observer 
manually draws an outline of the CME. All the pixels within 
this ROI is considered to contribute towards the CME mass. 
This method calculates the total mass more accurately than 
the sector method but requires each frame to be considered 
individually, and is therefore more time-consuming. 

3. Graduated cylindrical shell model (CCS) method. 
This method uses the forward modeling technique developed 
by Thernisien et al. [2009] and Thernisien [2011] to define 
the outer boundary of the CME. The different parameters in 
the model apart from the radial distance are manually cho- 
sen from comparing the model shape with both STEREO 
COR-2 images simultaneously. These parameters are then 
fixed and the radial distance is varied over the sequence of 
images as the CME propagates. The outer edge which is 
traced out by the model is used in the same manner as the 
ROI method defined above. This new CCS method is used 
in our paper. 

The CME enters the COR-2A camera at 21.24 UT, March 
24 2011 and approaches a mass of ~ 3 x 10 12 fcg which is of 
the order of magnitude of a typical CME [Vourlidas et al, 
2010]. The mass estimates as measured from the CCS model 
are displayed in figure 3 and show the typical increasing- 
mass time profile due to the CME entering the field of view. 
This method enables the CME mass calculations to be made 
in a semi-automated fashion similar to the sector method, 
while using a more reliable ROI method for tracing the 
CME contour. In our example the mass of the CME was 
slightly overestimated by < 5% compared to using method 
2. This is because we ensured that the entire CME struc- 
ture was enclosed. For our event, which may not always 
be the case, it was relatively simple as the CME had well 
defined boundaries. However, as the CME propagated the 
leading edge became slightly flattened in comparison to the 
idealised CCS model. This meant the nose of the CCS was 
progressively further into the heliosphere than the observa- 
tions. This meant a few extra pixel of mass were included 
but less than would have been encountered by method 1. It 
was found that the accuracy of the mass measurements com- 
pared to method 2 depended on ensuring that at least the 
entire area of the CME was included rather than minimising 
the surplus area of background solar wind. 

The CCS method used in this paper requires several pa- 
rameters to be defined for each image. The results of prop- 
agation direction and flux rope axis were compared to the 
fixed-phi method (from remote observations) and constant- 
alpha force-free flux rope model (from in situ data; see § 3.1 
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for more details), respectively. The results were consistent 
with each other. The independently estimated propagation 
directions were within 10° of each other and the three inde- 
pendently estimated flux rope axes were within 30° 

Another option for measuring the mass flow is to investi- 
gate a fixed location in the heliosphere by defining a small 
narrow rectangular box a few pixels wide (slit method). This 
method emulates a time series of mass over a fixed location 
such as a spacecraft at the LI position (see § 2.4). Cur- 
rently, the minimal image processing used in this study does 
not allow for an accurate estimate to be made at Lf . How- 
ever, future studies using more advanced image processing 
and better tempo-spatial resolution as expected from Solar 
Orbiter should provide the necessary data to advance the 
techniques in preparation for a possible mission to L5. 

2.3. Mass J-maps 

Originally developed for LASCO coronagraphic images, 
Sheeley et al. [1 999] calculated that a small plasma packet 
moving at uniform speed would have an apparent acceler- 
ation and deceleration which is dependant on 2 variables: 
the radial velocity, Vr, and the angle of propagation be- 
tween the CME-Sun-spacecraft, /3. Moreover, the observed 
acceleration profile (measured in time and elongation an- 
gle away from the Sun, a) is unique; therefore by using 
an optimisation routine [e.g. Savani et al, 2009, 20f0], an 
estimate of the propagation direction of the CME can be 
made. This technique (called fixed-phi method) has been 
shown to be much more useful over a large range of elon- 
gation angles [Williams et al, 2009; Davis et al, 20f0] and 
produced very effective results in tracking CMEs from the 
Sun to planetary systems where they were detected in situ 
[e.g. Davis et al., 2009; Rouillard et al., 2009; Mostl et al., 
2009]. As the premise on which this technique is built re- 
lies on the idea of a spatially narrow plasma packet and 
not a large 3-dimensional (3D) object travelling through the 
heliosphere, other attempts have been developed to miti- 
gate against some of the simplified assumptions [Kahler and 
Webb, 2007; Lugaz et al, 2009; Davies et al, 2012]. 

The propagation direction is estimated by measuring the 
elongation angle as a function of time. This is most con- 
veniently done by tracking a single feature within a CME 
through the field of view of all the cameras by using a time- 
elongation map (J- map). As a plasma packet can be as- 
sumed to propagate radially away from the Sun, the radial 
cuts used in J-maps can be varied to suite the CME direction 
if it is away from the ecliptic plane. In the case of our event 
we choose a P.A. of 97° from ST-A, which corresponds to 
the direction to Earth. Until this paper, these J-maps have 
been created by processing MSB images into a running dif- 
ference sequence, as shown in the top panel of Figure 4. 
The blue crosses are manually chosen from the image and 
are the points used in estimating the speed and propagation 
direction of the CME. Using the fixed-phi method the ra- 
dial speed and propagation direction was estimated as 371 
km/s and 80°, respectively. The bottom panel displays the 
same time-elongation tracks but created with a sequence of 
mass images. As can be seen between the two versions, the 
white light J-map is clearer at identifying the propagation 
of the CME. This is partly because further smoothing and 
image processing was undertaken on the white-light J-map. 
It is hoped that further studies into mass J-maps will be 
used to provide a time series estimate of the mass propagat- 
ing over the Earth in a process similar to the slit method 
described above. This is because the slit method would ef- 
fectively represent a horizontal line along a mass J-map. 
Although the tracks in the mass J-map displayed in our pa- 
per is 'noisy', future studies may implement advanced image 
processing techniques that are currently under development 
[e.g. Howard et al, 2012] to find interesting discoveries. 

The tracking of the mass estimates can be used with a 
minimal number of assumptions about the expansion pro- 
cess to estimate the density of the CME. With improvements 



to the image processing this information may be used as in- 
puts to space weather forecasting models, instead of using 
LI in situ data. Although this will clearly be less reliable 
than LI data itself, it has the big advantage of being mea- 
sured remotely and of the order of ~ 48 hours in advance. 
This may prove to be a significant improvement for our fore- 
casting capabilities. 

Figure 5 displays the mass J-map along with the in situ 
measurements of density and velocity at LI. The track, 
which was made from the white-light J-map, is overplotted 
onto the mass J-map along with a dashed line to show the 
position of Earth. The track clearly intercepts the position 
of Earth at the same time as the CME-driven shock (and 
the associated density increase from the sheath) arrives at 
LI. The shock arrival at LI displays a sudden increase in 
momentum flux which then compresses the Earth's magne- 
tosphere (see § 4). 

2.4. Mass time series 

Under the premise that the CME momentum can be 
tracked and possibly be used as early solar wind input into 
space weather simulations, it is important to observe how 
a typical mass time series may look like and how it com- 
pares to the in situ density profile currently being used as 
simulation inputs. Figure 5 shows a normalised times se- 
ries of the mass measurements (green dashed curve) taken 
at 3.8° elongation from ST-A, which is within the COR-2 
field of view. The mass measurements were taken using the 
slit method at a plane of sky distance between 13.7-14Rs for 
each frame. The data was then ballistically time-shifted to 
the LI position by assuming the CME travelled at 370km/s 
and was linearly expanding [Owens et al, 2005] so that the 
trailing edge of the mass measurements propagated 30km/s 
slower than the sheath leading edge (see also Figure 6 and 7). 
As our event initially propagates with a slow speed we are 
able to assume the CME was swept into the solar wind and 
merely advected out to 1 AU [Siscoe and Schwenn, 2006]. 
However for faster CMEs it is important to consider decel- 
eration due to drag effects that change the arrival times at 
Earth [Gopalswamy et al, 2000; Cargill, 2004]. The units 
of mass displayed have been normalised to suit the density 
profile. This has simply been carried out by dividing the 
measurements by 1.5 x 10 g.cm . A simple method of an 
expanding volume during the CME's propagation is used to 
justify this value (see Appendix A for more details). 

Currently only minimally processed images are being used 
and even with our simplified propagation technique [Owens 
and Cargill, 2004; Lugaz and Kintner, 2012, e.g], the quali- 
tative profile of the mass measurements made from remote 
observations is remarkably similar to that observed in situ. 
It is clear that CMEs may undergo interactions during prop- 
agation to Earth by either solar wind distortions [e.g. Lugaz 
et al, 2008; Savani et al, 2010], deflections [Lugaz et al, 
2012; Wood et al, 2012] or possible rotations [e.g. Shiota 
et al, 2005; Vourlidas et al, 2011; Nieves- Chinchilla et al, 
2012]. This could be the cause of the larger predicted mass 
measurements within the CME when compared to the in 
situ values. Or more simply that the CME expanded more 
during the propagation and is more rarified by the time it 
reaches Earth. 

3. In situ observations 

The shock associated with our case study event was de- 
tected at LI at 15.09 March 29 2011 with a jump in velocity 
of ~ 70 km/s (the upstream speed was ~ 330 km/s). Figure 
6 shows the in situ paramaters of the solar wind during the 
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propagating interplanetary CME. The three vertical lines 
indicate the locations of the shock front (14.58, 29th), CME 
leading edge (23.39, 29th) and rear edge (09.52, 31st) re- 
spectively. The CME leading and rear edge were manually 
chosen by focusing on looking for a duration that includes 
a smooth rotation in magnetic field and a discontinuity on 
density. The focus on the smooth field rotation is in order to 
produce reliable results from an in situ flux rope fitting pro- 
cess. The magnetic field vectors are displayed in the RTN 
coordinate system such that the Normal (N) component is 
the measure of the magnetic field in the out-of ecliptic direc- 
tion (i.e. the North-South direction that is crucial in space 
weather predictions). This case study shows that the mag- 
netic field in the CME and sheath region is predominately 
northwardly directed, and strongly so in the earliest half of 
the CME just behind the leading edge. This CME topology 
was specially chosen for our analysis as it allows our study to 
isolate the space weather effects (e.g. potential strength of 
GICs) that are caused by the momentum flux and not from 
the resulting magnetic reconnection from a southwardly di- 
rected interplanetary magnetic field (IMF). 

3.1. Modelling results for March 2011 

Figure 6 displays the magnetic field profile of the opti- 
mised CAFF model (see Appendix B for details) within the 
top 6 panels as black curves. The top 3 represent the carte- 
sian vectors in RTN coordinate system and the other 3 rep- 
resent the field vectors in spherical coordinate system. The 
smooth rotation in the field is often more clearly seen in 
spherical coordinates while the importance of a southward 
Bz for predicting space weather events is better observed in 
cartesian coordinates. The orientation of the estimated flux 
rope axis direction is (0.4, —0.7, 0.6) in RTN and has a right 
handed chirality. The optimised parameters for the axial 
magnetic field and the impact parameter are 14.0 nT and 
0.1, respectively. This indicates the spacecraft travel close- 
to but slightly above the FR axis. The mean square error 
between the optimised model and the data, \, was < 0.1, 
which represents a good fit to the data [Lynch et at, 2003]. 
The model fits were also compared to a non-force-free ellip- 
tical model [Hidalgo and Nieves- Chinchilla, 2012] to show 
consistent results (e.g. the axis orientation were within 30° 
of each other). 

3.2. Dst index 

There is a variety of ground responses at Earth from space 
weather disturbances and measurements are often focused 
on different geographically localised processes/positions 
(e.g. due to local noon time, latitude or ground resistiv- 
ity). The disturbance time (Dst) index measures the hourly 
values of the horizontal component of the Earth's magnetic 
field averaged from four near-equatorial geomagnetic obser- 
vatories. The Dst index has historically been used as an ap- 
proximation of the global response to a space weather distur- 
bance. The fluctuations in the Dst closely relates to the ring 
current and the other current systems (including the mag- 
netotail current) within the terrestrial environment. The 
inverse proportionality relationship between the horizontal 
component of the magnetic field and the energy content of 
the ring current is known as the Dessler-Parker-Sckopke re- 
lation [Dessler and Parker, 1959; Sckopke, 1966]. 

The stereotypical time profile of the Dst index during a 
geomagnetic storm driven by a fast CME displays an initial 
positive sharp rise (called the sudden storm commencement, 
SSC) which defines the arrival of the leading shock front onto 
the magnetosphere, a drop in value to zero (called the storm 
onset, SO), and then the main phase is characterised by a 
large negative decrease which represents the period of strong 
southward Bz magnetic field. The sheath region between the 
leading edge of a CME and the shock front is often consid- 
ered to be related to the period between the SSC and the 



SO (initial phase) during a geomagnetic storm. However, 
the phases of a geomagnetic storm can be time shifted un- 
der conditions where a strong southward Bz persists within 
the sheath region. 

The case study CME analysed in this paper was cho- 
sen to better understand the significance of its momentum 
flux as a driver of geomagnetic storms. As such, we chose 
to investigate a CME with a predominately northward Bz 
field. This allowed us to isolate the observed geomagnetic 
disturbance and assume the disturbance is solely due to the 
momentum flux and not due to reconnection between the 
CME and magnetosphere. For this reason, the usual main 
phase of a storm as seen in the Dst is not seen in our event 
(Figure 7). However, Figure 7 shows a significant rise in 
the Dst value and sharp fall during the initial phase, which 
represents the location of a sudden increase in momentum 
flux from the CME sheath region. The sudden changes in 
the Dst shows that the momentum flux of the CME lead- 
ing edge is capable of making sudden changes to the ground 
magnetic field (i.e. cause a large dB/dt). 

4. Terrestrial response 

Sophisticated MHD simulations in 3D are becoming an in- 
creasingly effective tool for modelling solar wind transients 
such as CMEs and for predicting their geo-effectiveness at 
Earth. In this paper we have employed the Space Weather 
Modeling Framework [Toth et al, 2005, SWMF] package 
which was executed at the Community Coordinated Model- 
ing Center (CCMC) and operated at NASA Goddard Space 
Flight Center. The solar wind input data was chosen from 
the WIND spacecraft at LI and was ballistically mapped 
to the outer boundary of the BATSRUS magnetospheric 
MHD model that was coupled to the Rice Convection Model 
(RCM) in our simulations. The auroral conductances are 
driven by solar irradiance observations of F10.7 and field- 
aligned electric currents. 

Figure 8 displays a 2D cut of Earth's magnetosphere (a) 
before the arrival of the CME-associated shock front and 
(b) when the sheath region of the CME is travelling over 
Earth's bow shock. The colour table represents the den- 
sity of the plasma and a few selected magnetic field lines 
are drawn to help distinguish between the terrestrial and 
heliospheric systems. The vectors show the solar wind di- 
rection. In the 3 hours between the frames shown in Figure 
8 the bow shock is severely compressed from a location of 
approximately 16Rs to ~ lORs. This compression signifi- 
cantly increased the density within the magnetosheath and 
was due to the sudden arrival of a larger momentum flux 
from within the sheath region of the CME. 

Currently this study uses the in situ data measured at LI 
as the inputs into the SWMF simulations in order to predict 
realistic geomagnetic disturbances. However it is envisaged 
that with further development to the mass images, a rea- 
sonable density time-series calculated remotely may be used 
as an input into the BATSRUS code (see § 2.4). This may 
prove to be a valuable tool in improving any early warn- 
ing systems by being able to provide an observationally- 
predicted result that is at least ~ 24 hours earlier than is 
currently possible. 

4.1. Observations of magnetic fluctuations 

In order to better understand the geomagnetic effects 
driven by CMEs and in particular to predict the socio- 
economic impacts from GICs it is important to study the 
localised effects observed at specific locations on Earth dur- 
ing the arrival of the CME. In this paper we study the mo- 
mentum flux of CMEs and therefore choose to investigate 
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the equatorially based ground stations. The higher lati- 
tude stations are likely to display magnetic field fluctuations 
that can be partly associated with auroral magnetosphere- 
ionosphere dynamics, and are therefore not the focus of this 
paper. The shock arrival at LI (~ 15.00 UT) defines the 
start of the geomagnetic disturbance. We investigate two 
locations on Earth that are determined to be locally noon 
and midnight at the shock arrival. The Vassouras (VSS) 
station in Brazil is used as the locally noon station. The 
Kanoya (KNY) and Kakadu (KDU) stations in Japan and 
Australia, respectively, are studied as the locally midnight 
stations for our investigation. 

In order to investigate the possible effects that a CME 
may have on GICs, the time-derivative of the ground mag- 
netic field (dB/dt) must be studied. Figure 9 displays the 
dB/dt at both local times and shows significant fluctuations 
at the arrival of the shock onto Earth's system. The ver- 
tical dashed line displays the time the shock arrived at LI 
and therefore the delay of ~ 1 hour on the ground stations 
predominately represents the propagation time of the so- 
lar wind between LI and the bow shock. The data is dis- 
played with one minute temporal resolution in the cartesian 
geographic (GEO) coordinate system provided by INTER- 
MAGNET (www.intermagnet.org). 

4.2. Simulated dB/dt 

Quantifying and predicting ground magnetic field pertur- 
bations are vital to the space weather community. As such, 
the CCMC has developed a tool that is able to extract the 
ground magnetic field perturbations from the global MHD 
model outputs by integrating the results from the magne- 
tospheric and ionospheric current systems. In particular, a 
summation of four separate current systems is used to make 
the predictions: 1. the current system in the magnetosphere 
(and magnetotail) above 2.5 Re; 2. field aligned currents 
between 2.5 Re and 110 km; 3. Hall current from the iono- 
sphere; 4. Pederson currents in the ionosphere [Rastaet- 
ter et al, 2004; Pulkkinen et al, 2010b]. The performance 
and metric-based analyses between various modelling ap- 
proaches were part of the Geospace Environment Modeling 
(GEM) 2008 - 2009 challenge and is reported by Pulkkinen 
et al. [2011]. 

Figure 9 shows the predicted estimates of the ground field 
as purple curves over-plotted on the data. The simulated 
time series have been shifted earlier in time by 15 minutes 
which we attribute to a small uncertainty in the input data 
which resulted from ballistically shifting the input LI data 
to the outer boundary of the simulation. We show that 
the simulated data also displays a significant spike in the 
dB/dt at the storm commencement which coincides with 
a similar magnitude to the 'ground truth' observed by the 
magnetometer ground stations. Qualitatively, the biggest 
limitation in the simulated results appears to occur in the 
inaccuracy of the negative dBz/dt component seen during 
the time derivative 'spike' This could be due to the Bz com- 
ponent being especially sensitive to geomagnetic induction 
effects which are not taken into account in our simulation. 

It is worth noting that only mid and high latitude mag- 
netometer stations were included in the earlier GEM chal- 
lenges even though the Dst index addresses the low-latitude 
disturbances. This is due to the global MHD approach be- 
ing implemented in the first-principle models; they could 
only capture the ionospheric output at high latitudes by 
using the Biot-Savart law to integrate over the ionospheric 
electric currents system. In this paper we have alleviated 
such a constraint by coupling the global MHD models to 
the inner magnetospheric models and thereby capturing the 
ring current dynamics and magnetospheric current systems; 
therefore providing the required ionospheric response at low 
latitudes [Yu et al, 2010]. 

4.3. Geomagnetically induced currents, GICs 

Reliable estimates of GICs and the geoelectric field re- 
quires accurate knowledge of the local geological conditions 



as well as the dB/dt. As the global distribution of the con- 
ductivity from the surface to the upper mantle (depths of 
several hundred kilometers) is not well known, estimating 
GICs can not be arbitrarily made for any location on Earth. 
However, for local environments that have historically been 
susceptible to GIC events [e.g. Pulkkinen et al., 2005; Ngwira 
et al, 2008; Pulkkinen et al, 2010a; Porta et al, 2012], the 
ground structure is known and therefore this study could 
be replicated for larger geo-effective CME events and early 
warning predictions can be made as to their level of suscepti- 
bility. It is envisaged that as the ground structure for more 
locations around the globe become recognized, the frame- 
work presented in this paper will provide the necessary steps 
to improve GIC forecasting. Of course, a simple approach 
may be created to estimate the extrema of possible GICs 
by using realistic extreme ends of the conducting (British 
Columbia, Canada) and resistive (Quebec, Canada) ground 
structure [Pulkkinen et al, 2008]. But detailed investigation 
into these possibilities goes beyond the scope of this work. 

5. Discussion 

Currently observations of the solar magnetic fields are 
used as inputs for the background solar wind when simulat- 
ing the heliosphere for space weather predictions. However 
the CME itself is usually set to an approximate and generic 
size. Details of the CME that is included in the magne- 
tospheric simulations are only provided by measurements 
made in situ at LI, which is ~ 1 hour before its arrival 
at the Earth's bow shock. By using remote observations 
from coronagraphs (~ 15Rs), observational estimates of the 
CME can be made for both the Enlil heliospheric simulation 
and as early initial-attempt inputs to BATSRUS magneto- 
spheric simulations. The use of remote observations means 
that BATRUS simulations can be carried out ~ 48 hours 
prior to waiting for the CME to propagate to LI. 

The focus of this paper was to estimate the CME mass as 
a time series close to the Sun; which was then ballistically 
propagated to LI. While previous studies have been carried 
out to better estimate the arrival times of a CME [Owens 
and Cargill, 2004; Taktakishvili et al, 2009, 2011] and exten- 
sive plane-of-sky speed measurements have been made over 
a large number of CMEs [St. Cyr et al, 2000; Yashiro et al, 
2004]. These studies could be used to remotely estimate a 
velocity time series which can then be used to propagate 
the CME to LI. The ability to remotely estimate both the 
velocity and mass allows early prediction of the momentum 
flux arriving to Earth. 

Figure 5 and 7 display the remotely-observed mass time- 
series which have been artificially normalised to suit the in 
situ number density. The mass time-series were arbitrarily 
divided by 1.5 x 10 n </ in order to display a qualitative pro- 
file of the same order of magnitude. Clearly, if these mass 
time-series are to be used as inputs to early space weather 
simulations then an appropriate method to scale the mass 
estimates to a density value is required. Three basic meth- 
ods can be used: 1. Defining the CME's 3D volume by, 
for example, by using the GCS model [Thernisien et al, 
2006]. This volume can then be radially propagated out to 
1 AU where the density can be estimated. 2. Using the 
empirical formulas proposed by Vourlidas et al. [2010]. 3. 
Using a combination of white light images and off-limb spec- 
troscopy to directly measure the densities across the CME 
body. Future investigations may attempt to solve the most 
appropriate method with more events. Method 1 has been 
used in Appendix A to justify the order of magnitude used in 
this paper. With appropriate density and velocity estimates 
made remotely to first approximation, the results can be 
adjusted in an ensemble method for making space weather 
predictions along with variety of magnetic field estimates 
[e.g. Lin et al, 2000; Gopalswamy and Yashiro, 2011; Sa- 
vani et al, 2012b]. 
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6. Conclusions 

In this paper we show the first results for producing den- 
sity measurements from remote observations that are com- 
parable to their equivalent in situ time-series. Therefore 
this paper for the first time shows that a time series of data 
can be estimated remotely and be used to make forecasts of 
GICs at Earth. While a previous attempt by Pulkkinen et al. 
[2009] to forecast GICs from remote observations only used 
a generic pressure pulse for a CME within the Enlil model. 
We also confirm that the compression of the magnetosphere 
from the momentum flux of a CME is a significant variable 
in predicting geomagnetically induced currents. Our results 
are capable of producing qualitatively reliable estimates of 
the densities upstream of Earth's bow shock. These results 
indicate the possibility of using remote observations at a 
heliocentric distance of ~ 15Rs to estimate the solar wind 
density profile arriving at Earth and show that these esti- 
mates can be used as part of a preliminary early warning 
system for space weather predictions. 

This paper focuses on a case study event of a CME which 
displays a strong northward Bz component in the magnetic 
field rather than a fast geo-effective event. The CME prop- 
agated in the inner heliopshere between 25-31 March 2011 
and is observed in both ST-A and ST-B. The event allowed 
our study to attribute all the ground based effects at low 
latitudes to the compression of the magnetosphere and not 
to magnetic reconnection in the case of a southward directed 
Bz. The compression of the magnetosphere is found to be 
the result of the larger momentum flux (larger density and 
velocity) in the sheath region of the CME. 

Previous studies have estimated the mass of a CME as it 
travelled through the field of view by either manually mea- 
suring the region of interest (ROI) around the CME for each 
individual frame or by estimating it within a fixed sector and 
thereby assuming the inclusion of the additional solar wind 
mass is not significant in relation to the CME mass. Here, 
we improve the process by semi automating the ROI while 
refraining from the manual selections within each frame. We 
do this by using the graduated cylindrical shell model [Th- 
ernisien, 2011] to define the ROI and manually fixing the 
parameters by eye at the beginning and then by only varying 
the radial distance for each consecutive frame. 

To track the CME mass from the remote observations 
to LI we employed the J-map technique for a sequence of 
mass images in the STEREO data. The J-map technique 
has previously been extensively used with white light images 
to estimate the propagation direction and arrival time of a 
CME [e.g. Sheeley et al, 1999; Rouillard et al, 2009]. How- 
ever, we show that the mass J-map technique can be used 
to highlight the mass intensities travelling towards Earth. 
As mass images are estimated from base-difference frames 
rather than running differences, the background star field 
becomes significant in the HI- 2 field of view. In this paper 
we have chosen to use a minimal amount of image processing 
in order to estimate a base confidence level for the uncer- 
tainties in the mass values. Future studies may consider a 
comparison between the minimal processed images to more 
advanced techniques. 

As the mass images in the vicinity of Earth (HI-2 field of 
view) was not clearly visible in the images due to the inter- 
ference of the background star field, we chose to estimate the 
mass as a time series at a fixed heliocentric location within 
the COR-2A field of view (plane of sky distance of 14 Rs). 
The mass time series was then time shifted by assuming a 
leading edge speed of 400km/s and having a linear expan- 
sion profile such that the trailing edge of the mass profile 
was travelling at 370km/s. The qualitative comparison of 
our estimated mass profile to the in situ density measure- 
ments at LI are remarkably similar for the sheath region 
between the shock front and leading edge. Future studies 
to investigate CMEs may not necessarily have the capabili- 
ties to track the CME along the entire Sun-Earth line. As 
such, our ballistic propagation approach for the mass mea- 
surements from ~ 15Rs to LI may serve as a useful tool for 
estimating a time series for the mass. 



Appendix A: Density Normalisation 

In this paper we have chosen to artificially divide the 
time series of the mass measurements by a constant value 
of 1.5 x 10 11 ;;. This was used to convert a mass times series 
that was estimated remotely into a density time series which 
was later compared with the in situ measurements. Here we 
carry out some preliminary work to justify the number used 
with a simple volume expansion method. The work below 
is intended to provide an order of magnitude justification, 
however further work would benefit from a more detailed 
approach as suggested in § 5. 

We assume that a small volume (defined by the size of a 
few pixels in the COR 2A camera) measured at ~ lORs has 
a shape of a cylinder, such that the circular cross section is 
within the plane of sky and the length (Wo) is defined along 
the line of sight. First we assume that the length of the 
cylinder can be estimated from the CCS model (see Figure 
10 for details). At 15.24 UT on March 25 2011, the height 
of the legs (h) and the half angular width (a) from CCS 
as defined by Thernisien et al. [2006] was estimated as 9.78 
Rs and 32.5°, respectively. With trigonometry, we calcu- 
late that the cylindrical length should be 10.51 Rs. For this 
time, we noted that the plane of sky position of the CME 
nose was 10.5 Rs. The same process was carried out for 
an image at 17.24 UT to deduce that the heliocentric dis- 
tance of the CME was to first approximation equal to the 
cylindrical length, Wo- Therefore, at LI where the in situ 
measurements are made, we estimate the cylindrical length 
as Wi = 2WRs. 

To estimate the cross sectional area we assume the cross- 
section expands at a uniform speed of V ex = 30fcm/s (which 
we estimated from the in situ measurements). The total 
time of expansion, t, is the same as the propagation time and 
dependant on the bulk flow speed which we assume to be 
370 km/s. Assuming the propagation distance, R =200Rs, 
we find that the radius of the circular cross section at LI is 
16.22Rs and follows 



n = 



Vbuh 



x R 



(Al) 



Therefore our final volume element (Y f) can be estimated 



by, 



(A2) 



to give 5.85 x 10 37 cm 3 . 

It therefore follows that a first approximation for a nor- 
malisation constant (/to) is, 



K() = 



From the constants used in our example, Ko 



(A3) 



1.02 x 



io- n ff - 



cm . This is similar to the normalisation value 
of 1.5 x 10 11 used in this paper and is certainly of the correct 
order of magnitude. This small discrepancy might be solved 
with more advanced calculations. 



Appendix B: Flux rope fitting 
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Simple determination of parameters from a model fitting 
procedure is one of the best ways to quickly estimate the 
global properties of an interplanetary CME from in situ mea- 
surements. The first model to be successfully optimised to 
MCs was a constant a, force-free (CAFF) flux rope model 
by Lepping et al. [1990]. Since then, several other attempts 
have been made to improve the results between observa- 
tions and models [Owens et al, 2006; Mulligan and Russell, 
2001; Hidalgo et al, 2002; Hidalgo and Nieves-Chinchilla, 
2012; Owens et al, 2012]. We use a simplified modification 
of the Lepping et al. [1990] model. The magnetic field of 
the model is described by the cartesian form of the Bessel 
function for a force-free flux rope. In fitting the flux rope 
model to the data, the first step is to estimate the orienta- 
tion of the cylindrical axis. Then by determining the chiral- 
ity and varying both the distance of closest approach and 
axial magnetic field strength as free parameters within a 
computational code, we are able to make predictions of the 
ICME parameters. 

The method employed to determine the MC axis direc- 
tion involves a technique called minimum variance analysis, 
MVA [Sonnerup and Cahill, 1967; Steed et al., 2008]. MVA 
was originally developed to determine the local normal for 
a tangential discontinuity (TD) in the magnetic field. This 
method calculates the orthogonal set of vectors in which the 
variance of the magnetic field in question is at a minimum 
(el), maximum (e3) and intermediate (e2). For each of these 
eigenvectors there is a corresponding eigenvalue, often rep- 
resented as Ai,A3 and A2 respectively. 

For the purposes of investigating TDs the minimum eigen- 
vector calculated provides the normal to the plane. Here, the 
ratios of eigenvalues (Ai) and the other two are of concern as 
they indicate the reliability in the normal. For investigating 
the axial direction of a MC, we must concern ourselves with 
the intermediate direction, with MVA applied to the time 
period that contains the flux rope only. The axis of the MC 
therefore lies in the intermediate direction, but the distinc- 
tion of the field variance along the coordinate axes reduce 
as the spacecraft trajectory moves further from the central 
MC axis. 

When investigating the intermediate direction, the eigen- 
value ratios between A2/A1 and A3/A2 are of concern. Ideally 
both ratios should be > 10, but recent work using MVA have 
shown ratios of greater than 2 are adequate in the ICME 
context [Bothmer and Schwenn, 1998]. If the ratio of either 
A2/A1 or A3/A2 is small then the axial direction becomes 
ambiguous within the plane of the two eigenvectors. 

The scenario above is described for the ideal case, where 
the spacecraft passes through the centre of the flux rope. 
However, this is usually not true. The distance of closest 
approach, Yo, is an important parameter in determining the 
predicted field profile. It is defined in dimensionless units 
as the closest approach distance divided by the radius of 
the flux rope. Thus Yb varies between 1 and -1 for a space- 
craft travelling above and below the MC centre, respectively. 
Rees [2003] simulated idealised flux ropes with spacecraft 
trajectories away from the centre and then carried out MVA 
analysis over the time profiles. He found that the error in 
the axis angle increased with an increasing |Yo|. The results 
suggest that an impact parameter of > 0.5 should be dealt 
with caution as the error may be of the order of ~ 15°. The 
uncertainty in ascertaining the axis orientation [Al-Haddad 
et al., 2013] is one of the reasons why more complex models 
choose to begin their analysis with MVA and then later find 
the optimal direction by introducing 2 extra free variables. 

High frequency noise is a source of error for MVA in the 
interplanetary CME context. MVA analyses the variance 
in the magnetic field, therefore large noise fluctuations can 
yield incorrect eigenvectors. To limit this effect, 10 or 15 
minute averaged data are often used when fitting a mod- 
elled rope to the data. It has been shown for good events 
the time resolution between 10 and 60 minutes has a neg- 
ligible effect on the axis orientation [Lepping et al, 2003]. 



Subtly, the lower resolutions are often used for large events, 
with a reason to maintain the total number of data points 
within the MC to about N=40 [Lepping et al, 2006]. Al- 
though this is not essential, it is desirable when comparing 
various cases against each other in a statistical manner. 

Once the axis of the model MC is determined, the in situ 
magnetic field data are rotated into the frame of the MVA 
axes. A sequence of flux ropes are then created and com- 
pared to the data by analysing the mean square error (x) 
coefficient. The \ coefficient measures how well each rope 
fits to the data and is defined in the same manner as Lynch 
et al. [2003]. 

The time interval between the start and end of the MC 
is identified manually by eye. This defines the data range 
which is compared to the model and fixes the size of the 
MC model. By using the distance of closest approach, the 
axial field strength and chirality (a concept of rotation in 
the magnetic field, either clockwise or anticlockwise) as free 
parameters, various simulated MCs can be generated. Each 
simulated MC is compared to the data by measuring the \ 
coefficient. The simulated MC with the smallest x value is 
regarded as the best fit result. The free variables that cre- 
ate the model are varied by using a downhill simplex method 
developed by N elder and Mead [1965]; this is a non- linear 
optimisation routine designed to minimise the \ coefficients. 
This approach is faster than a standard grid search and 
can be easily manipulated to determine the free variables 
to larger significant figure. The optimised free parameters 
are the outputs of the model. These results are used again 
to re-create the optimised model. The field vectors from 
this model are then rotated into the spacecraft frame and 
plotted on top of the observed data for visual confirmation. 

It is worth noting here that the CAFF model described 
above assumes a static ICME. That is to say, the model 
plotted shows a time series obtained by taking a radial cut 
through a flux rope at a fixed time. Many ICMEs observed 
at terrestrial distances take ~ 24 hours to transit over a 
spacecraft, therefore it is important to note that the ICME 
seen at the end of the transit has evolved and is different 
from the same object at the beginning of the transit. A 
more accurate representation would be to include a time 
series of a flux rope evolving in time past a fixed point in 
space. 

Acknowledgments. NPS thanks Neil Shccley and Adam Sz- 
abo for their collaborative assistance. This research was also sup- 
ported by the NASA Living With a Star Jack Eddy Postdoctoral 
Fellowship Program, administered by the UCAR Visiting Scien- 
tist Programs and hosted by the Naval Research Laboratory. AV 
is supported by NASA contract S-136361-Y. The global MHD 
simulations used in this work were carried out at the Commu- 
nity Coordinated Modeling Center (CCMC) operated at NASA 
Goddard Space Flight Center. The authors wish to acknowl- 
edge Lutz Rastaetter, David Berrios and the rest of the CCMC 
staff for their generous support throughout the work discussed 
in the paper. The results presented in this paper rely on the 
data collected at Vassouras, Kakadu and Kanoya. We thank Ob- 
servatorio Nacional, Geoscience Australia and Japan Meterolog- 
ical Agency, for supporting its operation and INTERMAGNET 
for promoting high standards of magnetic observatory practice 
(www. intermagnct .org) . 



References 

Al-Haddad, N., T. Nieves-Chinchilla, N. P. Savani, C. Mostl, 
K. Marubashi, M. Hidalgo, I. I. Roussev, S. Poedts, and 
C. J. Farrugia (2013), Magnetic Field Configuration Mod- 
els and Reconstruction Methods for Interplanetary Coro- 
nal Mass Ejections, Solar Physics, pp. In Press,, doi: 
10.1029/JA090iA07p06387. 



SAVANI ET AL.: MOMENTUM FLUX OF CMES 



X- 9 



Anderson, C. W. I., L. J. Lanzerotti, and C. G. MacLennan 
(1974), Outage of the L4 system and the geomagnetic distur- 
bances of 4 August 1972, The Bell System Technical Journal, 
9, 1817. 

Bargatze, L. F., D. N. Baker, E. W. Hones, Jr., and R. L. McPher- 
ron (1985), Magnetospheric impulse response for many levels 
of geomagnetic activity, J. Geophys. Res., , 90, 6387-6394, 
doi: 10. 1029/ JA090iA07p06387. 

Billings, D. E. (1966), A guide to the solar corona. 

Borovsky, J. E, and M. H. Denton (2006), Differences be- 
tween CME-driven storms and CIR-driven storms, Journal 
of Geophysical Research (Space Physics), 111, A07S08, doi: 
10.1029/2005JA011447. 

Borovsky, J. E., M. F. Thomsen, and R. C. Elphic (1998), The 
driving of the plasma sheet by the solar wind, J. Geophys. 
Res., , 103, 17,617-17,640, doi:10.1029/97JA02986. 

Bothmcr, V., and R. Schwenn (1998), The structure and origin of 
magnetic clouds in the solar wind, Annales Geophysicae, 16, 
1-24, doi:10.1007/s00585-997-0001-x. 

Burlaga, L. F. (1988), Magnetic clouds and force-free fields 
with constant alpha, J. Geophys. Res., , 93, 7217-7224, doi: 
10.1029/JA093iA07p07217. 

Cargill, P. J. (2004), On the Aerodynamic Drag Force Acting 
on Interplanetary Coronal Mass Ejections, Solar Physics, 221, 
135-149, doi: 10. 1023/B:SOLA. 0000033366. 10725.a2. 

Committee On The Societal, and Economic Impacts Of Severe 
Space Weather Events (2008), Severe Space Weather Events- 
Understanding Societal and Economic Impacts: A Workshop 
Report, Tech. rep. 

Davies, J. A., R. A. Harrison, A. P. Rouillard, N. R. Sheeley, C. H. 
Perry, D. Bewsher, C. J. Davis, C. J. Eyles, S. R. Crothcrs, 
and D. S. Brown (2009), A synoptic view of solar transient 
evolution in the inner heliosphere using the Heliosphcric Im- 
agers on STEREO, Geophys. Res. Lett., , 36, L02102, doi: 
10. 1029/2008GL036182. 

Davies, J. A., R. A. Harrison, C. H. Perry, C. Mostl, N. Lugaz, 
T. Rollett, C. J. Davis, S. R. Crothers, M. Temmer, C. J. Eyles, 
and N. P. Savani (2012), A Self-similar Expansion Model for 
Use in Solar Wind Transient Propagation Studies, Astrophys. 
J., , 750, 23, doi:10.1088/0004-637X/750/l/23. 

Davis, C. J., J. A. Davies, M. Lockwood, A. P. Rouillard, C. J. 
Eyles, and R. A. Harrison (2009), Stereoscopic imaging of an 
Earth-impacting solar coronal mass ejection: A major mile- 
stone for the STEREO mission, Geophys. Res. Lett., , 36, 
L08102, doi: 10. 1029/2009GL038021 . 

Davis, C. J., J. Kennedy, and J. A. Davies (2010), Assessing 
the Accuracy of CME Speed and Trajectory Estimates from 
STEREO Observations Through a Comparison of Independent 
Methods, Solar Physics, 263, 209-222, doi:10.1007/sll207- 
010-9535-2. 

Denton, M. H., J. E. Borovsky, R. M. Skoug, M. F. Thom- 
son, B. Lavraud, M. G. Henderson, R. L. McPherron, J. C. 
Zhang, and M. W. Liemohn (2006), Geomagnetic storms 
driven by ICME- and CIR-dominated solar wind, Journal 
of Geophysical Research (Space Physics), 111, A07S07, doi: 
10.1029/2005JA011436. 

Dossier, A. J., and E. N. Parker (1959), Hydromagnctic theory of 
geomagnetic storms, J. Geophys. Res., , 64, 2239-2252, doi: 
10.1029/JZ064i012p02239. 

Eyles, C. J., R. A. Harrison, C. J. Davis, N. R. Waltham, B. M. 
Shaughnessy, H. C. A. Mapson-Menard, D. Bewsher, S. R. 
Crothcrs, J. A. Davies, G. M. Simnett, R. A. Howard, J. D. 
Moses, J. S. Newmark, D. G. Socker, J. -P. Halain, J.-M. De- 
fise, E. Mazy, and P. Rochus (2009), The Heliosphcric Imagers 
Onboard the STEREO Mission, Solar Physics, 254, 387-445, 
doi:10.1007/sll207-008-9299-0. 

Fukushima, N., and Y. Kamide (1973), Partial ring cur- 
rent models for worldwide geomagnetic disturbances., Re- 
views of Geophysics and Space Physics, 11, 795-853, doi: 
10.1029/RG011i004p00795. 

Gonzalez, W. D., and F. S. Mozer (1974), A quantitative model 
for the potential resulting from reconncction with an arbitrary 
interplanetary magnetic field, J. Geophys. Res., , 79, 4186- 
4194, doi: 10. 1029/ JA079i028p04186. 

Gopalswamy, N., and S. Yashiro (2011), The Strength and Ra- 
dial Profile of the Coronal Magnetic Field from the Standoff 
Distance of a Coronal Mass Ejection-driven Shock, Astrophys. 
J., , 736, L17, doi:10.1088/2041-8205/736/l/L17. 



Gopalswamy, N., A. Lara, R. P. Lopping, M. L. Kaiser, 
D. Berdichevsky, and O. C. St. Cyr (2000), Interplanetary 
acceleration of coronal mass ejections, Geophys. Res. Lett., , 
27, 145-148, doi:10.1029/1999GL003639. 

Hayes, A. P., A. Vourlidas, and R. A. Howard (2001), Deriving 
the Electron Density of the Solar Corona from the Inversion 
of Total Brightness Measurements, Astrophys. J., , 548, 1081— 
1086, doi:10.1086/319029. 

Hidalgo, M. A., and T. Nieves-Chinchilla (2012), A Global Mag- 
netic Topology Model for Magnetic Clouds. I., Astrophys. J., 
, 748, 109, doi:10.1088/0004-637X/748/2/109. 

Hidalgo, M. A., T. Nieves-Chinchilla, and C. Cid (2002), El- 
liptical cross-section model for the magnetic topology of 
magnetic clouds, Geophys. Res. Lett., , 29(13), 1637, doi: 
10.1029/2001GL013875. 

Hildner, E., J. T. Gosling, R. M. MacQueen, R. H. Munro, A. I. 
Poland, and C. L. Ross (1975), The large coronal transient of 
10 June 1973. I - Observational description, Solar Physics, 4%, 
163-177, doi:10.1007/BF00153293. 

Howard, R. A., N. R. Sheeley, Jr., D. J. Michels, and M. J. 
Koomen (1985), Coronal mass ejections - 1979-1981, J. Geo- 
phys. Res., , 90, 8173-8191, doi:10.1029/JA090iA09p08173. 

Howard, R. A., J. D. Moses, A. Vourlidas, J. S. Newmark, D. G. 
Socker, S. P. Plunkctt, C. M. Korendyke, J. W. Cook, A. Hur- 
ley, J. M. Davila, W. T. Thompson, O. C. St Cyr, E. Mentzell, 
K. Mchalick, J. R. Lemon, J. P. Wuelser, D. W. Duncan, 
T. D. Tarbell, C. J. Wolfson, A. Moore, R. A. Harrison, N. R. 
Waltham, J. Lang, C. J. Davis, C. J. Eyles, H. Mapson- 
Menard, G. M. Simnett, J. P. Halain, J. M. Dcfise, E. Mazy, 
P. Rochus, R. Mercier, M. F. Ravet, F. Delmotte, F. Auchere, 
J. P. Delaboudinicrc, V. Bothmer, W. Dcutsch, D. Wang, 
N. Rich, S. Cooper, V. Stephens, G. Maahs, R. Baugh, D. Mc- 
Mullin, and T. Carter (2008), Sun Earth Connection Coronal 
and Heliosphcric Investigation (SECCHI), Space Science Re- 
view, 136, 67-115, doi:10.1007/sll214-008-9341-4. 

Howard, T. A., and C. E. DeForest (2012), Inner Heliosphcric 
Flux Rope Evolution via Imaging of Coronal Mass Ejections, 
Astrophys. J., , 746, 64, doi:10.1088/0004-637X/746/l/64. 

Howard, T. A., C. E. DeForest, and A. A. Reinard (2012), White- 
light Observations of Solar Wind Transients and Compari- 
son with Auxiliary Data Sets, Astrophys. J., , 754, 102, doi: 
10. 1088 /0004-637X/754/2/102. 

Huttuncn, K. E. J., S. P. Kilpua, A. Pulkkincn, A. Viljanen, 
and E. Tanskanen (2008), Solar wind drivers of large geo- 
magnetically induced currents during the solar cycle 23, Space 
Weather, 6, S10002, doi:10.1029/2007SW000374. 

Jones, G. H., A. Rccs, A. Balogh, and R. J. Forsyth (2002), 
The draping of heliosphcric magnetic fields upstream of coro- 
nal mass ejecta, Geophys. Res. Lett., , 29(11), 1520, doi: 
10.1029/2001GL014110. 

Kahler, S. W., and D. F. Webb (2007), V arc interplanetary coro- 
nal mass ejections observed with the Solar Mass Ejection Im- 
ager, Journal of Geophysical Research (Space Physics), 112, 
A09103, doi: 10. 1029/2007JA012358. 

Kaiser, M. L., T. A. Kucera, J. M. Davila, O. C. St. Cyr, 
M. Guhathakurta, and E. Christian (2008), The STEREO 
Mission: An Introduction, Space Science Review, 136, 5—16, 
doi:10.1007/sll214-007-9277-0. 

Kataoka, R., and A. Pulkkinen (2008), Geomagnctically in- 
duced currents during intense storms driven by coronal 
mass ejections and corotating interacting regions, Journal 
of Geophysical Research (Space Physics), 113, A03S12, doi: 
10.1029/2007JA012487. 

Kataoka, R., S. Watari, N. Shimada, H. Shimazu, and 
K. Marubashi (2005), Downstream structures of interplanetary 
fast shocks associated with coronal mass ejections, Geophys. 
Res. Lett., , 32, L12103, doi:10.1029/2005GL022777. 

Lanzerotti, L. J. (1992), Comment on 'Great magnetic storms' 
by Tsurutani ct al, Geophys. Res. Lett., , 19, 1991, doi: 
10.1029/92GL02238. 

Lavraud, B., and J. E. Borovsky (2008), Altered solar wind- 
magnetosphere interaction at low Mach numbers: Coro- 
nal mass ejections, Journal of Geophysical Research (Space 
Physics), 113, A00B08, doi:10.1029/2008JA013192. 

Lavraud, B., and V. K. Jordanova (2007), Modeling the effects of 
cold-dense and hot-tenuous plasma sheet on proton ring cur- 
rent energy and peak location, Geophys. Res. Lett., , 34, 2102, 
doi: 10. 1029/2006GL027566. 



X - 10 



SAVANI ET AL.: MOMENTUM FLUX OF CMES 



Lavraud, B., M. F. Thomson, J. E. Borovsky, M. H. Denton, and 
T. I. Pulkkincn (2006), Magnetosphere preconditioning under 
northward IMF: Evidence from the study of coronal mass ejec- 
tion and corotating interaction region geoeffectiveness, Jour- 
nal of Geophysical Research (Space Physics), 111, A09208, 
doi:10.1029/2005JA011566. 

Lepping, R. P., L. F. Burlaga, and J. A. Jones (1990), 
Magnetic field structure of interplanetary magnetic clouds 
at 1 AU, J. Geophys. Res., , 95, 11,957-11,965, doi: 
10.1029/ JA095iA08pll957. 

Lepping, R. P., D. B. Berdichevsky, and T. J. Ferguson (2003), 
Estimated errors in magnetic cloud model fit parameters 
with force- free cylindrically symmetric assumptions, Jour- 
nal of Geophysical Research (Space Physics), 108, 1356, doi: 
10.1029/2002JA009657. 

Lepping, R. P., D. B. Berdichevsky, C.-C. Wu, A. Szabo, 
T. Narock, F. Mariani, A. J. Lazarus, and A. J. Quivers (2006), 
A summary of WIND magnetic clouds for years 1995-2003: 
model-fitted parameters, associated errors and classifications, 
Annates Geophysicae, 24, 215-245, doi:10.5194/angeo- 24-215- 
2006. 

Liemohn, M. W., J. U. Kozyra, M. F. Thomsen, J. L. Roeder, 
G. Lu, J. E. Borovsky, and T. E. Cayton (2001), Domi- 
nant role of the asymmetric ring current in producing the 
stormtime Dst*, J. Geophys. Res., , 106, 10,883-10,904, doi: 
10.1029/2000JA000326. 

Lin, H., M. J. Penn, and S. Tomczyk (2000), A New Precise Mea- 
surement of the Coronal Magnetic Field Strength, Astrophys. 
J., , 541, L83-L86, doi:10.1086/312900. 

Liu, C.-M., L.-G. Liu, R. Pirjola, and Z.-Z. Wang (2009), 
Calculation of gcomagnctically induced currents in mid- to 
low-latitude power grids based on the plane wave method: 
A preliminary case study, Space Weather, 7, S04005, doi: 
10.1029/2008SW000439. 

Lopez, R. E., M. Wiltberger, S. Hernandez, and J. G. Lyon 
(2004), Solar wind density control of energy transfer to 
the magnctosphere, Geophys. Res. Lett., , 31, L08804, doi: 
10.1029/2003GL018780. 

Lugaz, N., and P. Kintner (2012), Effect of Solar Wind Drag on 
the Determination of the Properties of Coronal Mass Ejec- 
tions from Heliospheric Images, Solar Physics, p. 47, doi: 
10.1007/sll207-012-9948-l. 

Lugaz, N., W. B. Manchester, IV, I. I. Roussev, and T. I. Gom- 
bosi (2008), Observational evidence of CMEs interacting in the 
inner hcliosphcrc as inferred from MHD simulations, Journal 
of Atmospheric and Solar-Terrestrial Physics, 70, 598-604, 
doi:10.1016/j.jastp.2007.08.033. 

Lugaz, N., A. Vourlidas, and I. I. Roussev (2009), Deriving the 
radial distances of wide coronal mass ejections from elonga- 
tion measurements in the heliosphcrc - application to CME- 
CME interaction, Annales Geophysicae, 27, 3479-3488, doi: 
10.5194/angeo-27-3479-2009. 

Lugaz, N., C. J. Farrugia, J. A. Davies, C. Mostl, C. J. Davis, 

I. I. Roussev, and M. Temmer (2012), The Deflection of the 
Two Interacting Coronal Mass Ejections of 2010 May 23-24 
as Revealed by Combined in Situ Measurements and Helio- 
spheric Imaging, Astrophys. J., , 759, 68, doi:10. 1088/0004- 
637X/759/1/68. 

Lynch, B. J., T. H. Zurbuchen, L. A. Fisk, and S. K. Antio- 
chos (2003), Internal structure of magnetic clouds: Plasma 
and composition, Journal of Geophysical Research (Space 
Physics), 108, 1239, doi:10.1029/2002JA009591. 

Miyoshi, Y., and R. Kataoka (2005), Ring current ions and radia- 
tion belt electrons during geomagnetic storms driven by coro- 
nal mass ejections and corotating interaction regions, Geophys. 
Res. Lett., , 32, L21105, doi:10.1029/2005GL024590. 

Mostl, C, C. J. Farrugia, M. Temmer, C. Miklenic, A. M. Veronig, 
A. B. Galvin, M. Leitner, and H. K. Biernat (2009), Linking 
Remote Imagery of a Coronal Mass Ejection to Its In Situ 
Signatures at 1 AU, Astrophys. J., , 705, L180-L185, doi: 
10.1088/0004-637X/705/2/L180. 

Mulligan, T., and C. T. Russell (2001), Multispacecraft mod- 
eling of the flux rope structure of interplanetary coronal 
mass ejections: Cylindrically symmetric versus nonsymmet- 
ric topologies, J. Geophys. Res., , 106, 10,581-10,596, doi: 
10.1029/2000JA900170. 

Nakagawa, T., A. Nishida, and T. Saito (1989), Planar magnetic 
structures in the solar wind, J. Geophys. Res., , 94, 11,761- 

II, 775, doi:10.1029/JA094iA09pll761. 



Nelder, J. A., and R. Mead (1965), A simplex method for 
function minimization, Computer Journal, 7, 308-313, doi: 
10.1093/comjnl/7.4.308. 

Ngwira, C. M., A. Pulkkincn, L.-A. McKinncll, and P. J. Cilliers 

(2008) , Improved modeling of geomagnetically induced cur- 
rents in the South African power network, Space Weather, 6, 
S11004, doi:10.1029/2008SW000408. 

Nieves-Chinchilla, T., R. Colaninno, A. Vourlidas, A. Szabo, 
R. P. Lepping, S. A. Boardscn, B. J. Anderson, and H. Korth 
(2012), Remote and in situ observations of an unusual Earth- 
directed coronal mass ejection from multiple viewpoints, Jour- 
nal of Geophysical Research (Space Physics), 117, A06106, 
doi:10. 1029/2011 JA017243. 

O'Brien, T. P., and R. L. McPherron (2000), Evidence against 
an independent solar wind density driver of the terrestrial 
ring current, Geophys. Res. Lett., , 27, 3797-3800, doi: 
10.1029/2000GL012125. 

Owens, M., and P. Cargill (2004), Predictions of the arrival time 
of Coronal Mass Ejections at 1AU: an analysis of the causes of 
errors, Annales Geophysicae, 22, 661-671, doi:10.5194/angco- 
22-661-2004. 

Owens, M. J., P. J. Cargill, C. Pagel, G. L. Siscoe, and N. U. 
Crooker (2005), Characteristic magnetic field and speed prop- 
erties of interplanetary coronal mass ejections and their sheath 
regions, Journal of Geophysical Research (Space Physics), 110, 
A01105, doi:10.1029/2004JA010814. 

Owens, M. J., V. G. Merkin, and P. Riley (2006), A kinemati- 
cally distorted flux rope model for magnetic clouds, Journal 
of Geophysical Research (Space Physics), 111, A03104, doi: 
10.1029/2005JA011460. 

Owens, M. J., P. Demoulin, N. P. Savani, B. Lavraud, and A. Ruf- 
fenach (2012), Implications of Non-cylindrical Flux Ropes for 
Magnetic Cloud Reconstruction Techniques and the Interpre- 
tation of Double Flux Rope Events, Solar Physics, 278, 435- 
446, doi:10.1007/sll207-012-9939-2. 

Pirjola, R. (2002), Review On The Calculation Of Surface Elec- 
tric And Magnetic Fields And Of Geomagnetically Induced 
Currents In Ground-Based Technological Systems, Surveys in 
Geophysics, 23, 71-90, doi:10.1023/A:1014816009303. 

Poland, A. I., R. A. Howard, M. J. Koomen, D. J. Michels, and 
N. R. Sheeley, Jr. (1981), Coronal transients near sunspot max- 
imum, Solar Physics, 69, 169-175, doi:10.1007/BF00151264. 

Pulkkincn, A., S. Lindahl, A. Viljanen, and R. Pirjola (2005), 
Geomagnetic storm of 29-31 October 2003: Geomagnetically 
induced currents and their relation to problems in the Swedish 
high-voltage power transmission system, Space Weather, 3, 
S08C03, doi:10.1029/2004SW000123. 

Pulkkincn, A., A. Viljanen, and R. Pirjola (2006), Estimation 
of geomagnetically induced current levels from different input 
data, Space Weather, 4, S08005, doi:10.1029/2006SW000229. 

Pulkkincn, A., R. Pirjola, and A. Viljanen (2007), Determina- 
tion of ground conductivity and system parameters for optimal 
modeling of geomagnetically induced current flow in techno- 
logical systems, Earth, Planets, and Space, 59, 999-1006. 

Pulkkinen, A., R. Pirjola, and A. Viljanen (2008), Statistics of ex- 
treme geomagnetically induced current events, Space Weather, 

6, S07001, doi:10.1029/2008SW000388. 

Pulkkincn, A., A. Taktakishvili, D. Odstrcil, and W. Jacobs 

(2009) , Novel approach to geomagnetically induced current 
forecasts based on remote solar observations, Space Weather, 

7, S08005, doi:10.1029/2008SW000447. 

Pulkkincn, A., R. Kataoka, S. Watari, and M. Ichiki (2010a), 
Modeling geomagnetically induced currents in Hokkaido, 
Japan, Advances in Space Research, 4®, 1087-1093, doi: 
10.1016/j.asr.2010.05.024. 

Pulkkincn, A., L. Rastatter, M. Kuznctsova, M. Hesse, A. Ri- 
dley, J. Raeder, H. J. Singer, and A. Chulaki (2010b), Sys- 
tematic evaluation of ground and geostationary magnetic field 
predictions generated by global magnetohydrodynamic mod- 
els, Journal of Geophysical Research (Space Physics), 115, 
A03206, doi:10.1029/2009JA014537. 

Pulkkincn, A., M. Kuznetsova, A. Ridley, J. Raeder, A. Va- 
pirev, D. Weimer, R. S. Weigel, M. Wiltberger, G. Mill- 
ward, L. Rastatter, M. Hesse, H. J. Singer, and A. Chu- 
laki (2011), Geospace Environment Modeling 2008-2009 Chal- 
lenge: Ground magnetic field perturbations, Space Weather, 
9, S02004, doi:10.1029/2010SW000600. 



SAVANI ET AL.: MOMENTUM FLUX OF CMES 



X - 11 



Rastactter, L., M. M. Kuznctsova, M. Hesse, D. Dezeeuw, A. Rid- 
ley, T. Gombosi, J. Dorelli, and J. Raedcr (2004), Energy Flow 
from the Solar Wind Through Magnetosphere and Ionosphere 
in Global MHD Models, AGU Spring Meeting Abstracts, p. A5. 

Rouillard, A. P., J. A. Davies, R. J. Forsyth, A. Rees, C. J. 
Davis, R. A. Harrison, M. Lockwood, D. Bewsher, S. R. 
Crothers, C. J. Eyles, M. Hapgood, and C. H. Perry (2008), 
First imaging of corotating interaction regions using the 
STEREO spacecraft, Geophys. Res. Lett., , 35, L10110, doi: 
10.1029/2008GL033767. 

Rouillard, A. P., J. A. Davies, R. J. Forsyth, N. P. Savani, N. R. 
Sheeley, A. Thernisien, T.-L. Zhang, R. A. Howard, B. An- 
derson, C. M. Carr, S. Tsang, M. Lockwood, C. J. Davis, 
R. A. Harrison, D. Bewsher, M. Franz, S. R. Crothers, C. J. 
Eyles, D. S. Brown, I. Whittaker, M. Hapgood, A. J. Coates, 
G. H. Jones, M. Grande, R. A. Frahm, and J. D. Winning- 
ham (2009), A solar storm observed from the Sun to Venus 
using the STEREO, Venus Express, and MESSENGER space- 
craft, Journal of Geophysical Research (Space Physics), 114, 
A07106, doi:10.1029/2008JA014034. 

Savani, N. P., A. P. Rouillard, J. A. Davies, M. J. Owens, R. J. 
Forsyth, C. J. Davis, and R. A. Harrison (2009), The radial 
width of a Coronal Mass Ejection between 0.1 and 0.4 AU es- 
timated from the Hcliosphcric Imager on STEREO, Annates 
Geophysicae, 27, 4349-4358, doi:10.5194/angeo-27-4349-2009. 

Savani, N. P., M. J. Owens, A. P. Rouillard, R. J. Forsyth, and 
J. A. Davies (2010), Observational Evidence of a Coronal Mass 
Ejection Distortion Directly Attributable to a Structured So- 
lar Wind, Astrophys. J., , 714, L128-L132, doi:10.1088/2041- 
8205/714/1/L128. 

Savani, N. P., M. J. Owens, A. P. Rouillard, R. J. Forsyth, K. Ku- 
sano, D. Shiota, R. Kataoka, L. Jian, and V. Bothmcr (2011), 
Evolution of Coronal Mass Ejection Morphology with Increas- 
ing Heliocentric Distance. II. In Situ Observations, Astrophys. 
J., , 732, 117, doi:10.1088/0004-637X/732/2/117. 

Savani, N. P., J. A. Davies, C. J. Davis, D. Shiota, A. P. Rouillard, 
M. J. Owens, K. Kusano, V. Bothmer, S. P. Bamford, C. J. 
Lintott, and A. Smith (2012a), Observational Tracking of the 
2D Structure of Coronal Mass Ejections Between the Sun and 
1 AU, Solar Physics, 279, 517-535, doi:10.1007/sll207-012- 
0041-6. 

Savani, N. P., D. Shiota, K. Kusano, A. Vourlidas, and N. Lugaz 
(2012b), A Study of the Heliocentric Dependence of Shock 
Standoff Distance and Geometry using 2.5D Magnetohydro- 
dynamic Simulations of Coronal Mass Ejection Driven Shocks, 
Astrophys. J., , 759, 103, doi:10.1088/0004-637X/759/2/103. 

Sckopke, N. (1966), A general relation between the energy of 
trapped particles and the disturbance field near the Earth, J. 
Geophys. Res., , 71, 3125-3130, doi:10.1029/JZ071i013p03125. 

Sheeley, N. R., J. H. Walters, Y.-M. Wang, and R. A. Howard 

(1999) , Continuous tracking of coronal outflows: Two kinds of 
coronal mass ejections, J. Geophys. Res., , 104, 24,739-24,768, 
doi:10.1029/1999JA900308. 

Shiokawa, K., T. Ogawa, and Y. Kamide (2005), Low- 
latitude auroras observed in Japan: 1999-2004, Journal of 
Geophysical Research (Space Physics), 110, A05202, doi: 
10.1029/2004JA010706. 

Shiota, D., H. Isobc, P. F. Chen, T. T. Yamamoto, T. Sakajiri, 
and K. Shibata (2005), Self-Consistent Magnctohydrodynamic 
Modeling of a Coronal Mass Ejection, Coronal Dimming, and 
a Giant Cusp-shaped Arcade Formation, Astrophys. J., , 634, 
663-678, doi:10. 1086/496943. 

Siscoe, G., and R. Schwcnn (2006), CME Disturbance Forecast- 
ing, Space Science Review, 123, 453-470, doi:10.1007/sll214- 
006-9024-y. 

Siscoe, G. L., G. M. Erickson, B. U. 6. Sonncrup, N. C. May- 
nard, J. A. Schoendorf, K. D. Siebert, D. R. Weimer, W. W. 
White, and G. R. Wilson (2002), Hill model of transpolar po- 
tential saturation: Comparisons with MHD simulations, Jour- 
nal of Geophysical Research (Space Physics), 107, 1075, doi: 
10.1029/2001JA000109. 

Sonncrup, B. U. O., and L. J. Cahill, Jr. (1967), Magnetopausc 
Structure and Attitude from Explorer 12 Observations, J. Geo- 
phys. Res., , 72, 171, doi:10.1029/JZ072i001p00171. 

St. Cyr, O. C, S. P. Plunkctt, D. J. Michels, S. E. Paswaters, 
M. J. Koomen, G. M. Simnett, B. J. Thompson, J. B. Gur- 
man, R. Schwenn, D. F. Webb, E. Hildner, and P. L. Lamy 

(2000) , Properties of coronal mass ejections: SOHO LASCO 
observations from January 1996 to June 1998, J. Geophys. 
Res., , 105, 18,169-18,186, doi:10.1029/1999JA000381. 



Steed, K., C. J. Owen, L. K. Harra, L. M. Green, S. Dasso, A. P. 
Walsh, P. Demoulin, and L. van Driel-Gesztelyi (2008), Locat- 
ing the solar source of 13 April 2006 magnetic cloud, Annales 
Geophysicae, 26, 3159-3168, doi:10.5194/angeo-26-3159-2008. 

Stewart, R. T., M. K. McCabe, M. J. Koomen, R. T. Hansen, 
and G. A. Dulk (1974), Observations of Coronal Disturbances 
from 1 to 9 Rsun- I: First Event of 1973, January 11, Solar 
Physics, 36, 203-217, doi:10.1007/BF00151561. 

Subramanian, P., and A. Vourlidas (2007), Energetics of solar 
coronal mass ejections, Astron. Astrophys. , , 467, 685-693, 
doi:10. 1051/0004-6361:20066770. 

Taktakishvili, A., M. Kuznetsova, P. MacNeice, M. Hesse, 
L. Rastatter, A. Pulkkinen, A. Chulaki, and D. Odstrcil 
(2009), Validation of the coronal mass ejection predictions at 
the Earth orbit estimated by ENLIL heliosphere cone model, 
Space Weather, 7, S03004, doi:10.1029/2008SW000448. 

Taktakishvili, A., A. Pulkkinen, P. MacNeice, M. Kuznetsova, 
M. Hesse, and D. Odstrcil (2011), Modeling of coronal mass 
ejections that caused particularly large geomagnetic storms us- 
ing ENLIL heliosphere cone model, Space Weather, 9, S06002, 
doi:10.1029/2010SW000642. 

Thernisien, A. (2011), Implementation of the Graduated Cylin- 
drical Shell Model for the Three-dimensional Reconstruction 
of Coronal Mass Ejections, Astrophys. J. (Supp.), , 194, 33, 
doi: 10. 1088/0067-0049/194/2/33. 

Thernisien, A., A. Vourlidas, and R. A. Howard (2009), Forward 
Modeling of Coronal Mass Ejections Using STEREO/SECCHI 
Data, Solar Physics, 256, 111-130, doi:10.1007/sll207-009- 
9346-5. 

Thernisien, A. F. R., R. A. Howard, and A. Vourlidas (2006), 
Modeling of Flux Rope Coronal Mass Ejections, Astrophys. 
J., , 652, 763-773, doi: 10. 1086/508254. 

Thomsen, M. F., J. E. Borovsky, D. J. McComas, R. C. El- 
phic, and S. Maurice (1998), The magnctosphcric response to 
the CME passage of January 10-11, 1997, as seen at geosyn- 
chronous orbit, Geophys. Res. Lett., , 25, 2545-2548, doi: 
10.1029/98GL00514. 

Torta, J. M., L. Serrano, J. R. Regue, A. M. Sanchez, and 
E. Roldan (2012), Gcomagnetically induced currents in a 
power grid of northeastern Spain, Space Weather, 10, S06002, 
doi:10.1029/2012SW000793. 

Toth, G., I. V. Sokolov, T. I. Gombosi, D. R. Chesney, C. R. 
Clauer, D. L. De Zeeuw, K. C. Hansen, K. J. Kane, W. B. 
Manchester, R. C. Oehmke, K. G. Powell, A. J. Ridley, I. I. 
Roussev, Q. F. Stout, O. Volberg, R. A. Wolf, S. Sazykin, 
A. Chan, B. Yu, and J. Kota (2005), Space Weather Mod- 
cling Framework: A new tool for the space science commu- 
nity, Journal of Geophysical Research (Space Physics), 110, 
A12226, doi:10.1029/2005JA011126. 

Tsurutani, B. T., Y. T. Lee, W. D. Gonzalez, and F. Tang 
(1992a), Great magnetic storms, Geophys. Res. Lett., , 19, 
73-76, doi:10.1029/91GL02783. 

Tsurutani, B. T., W. D. Gonzalez, F. Tang, Y. T. Lee, M. Okada, 
and D. Park (1992b), Reply to L. J. Lanzcrotti: Solar wind 
RAM pressure corrections and an estimation of the efficiency 
of viscous interaction, Geophys. Res. Lett., , 19, 1993-1994, 
doi: 10. 1029/92GL02239. 

Viljancn, A., H. Ncvanlinna, K. Pajunpaa, and A. Pulkkinen 
(2001), Time derivative of the horizontal geomagnetic field 
as an activity indicator, Annales Geophysicae, 19, 1107-1118, 
doi:10.5194/angeo-19-1107-2001. 

Vourlidas, A., and R. A. Howard (2006), The Proper Treatment 
of Coronal Mass Ejection Brightness: A New Methodology and 
Implications for Observations, Astrophys. J., , 642, 1216-1221, 
doi:10.1086/501122. 

Vourlidas, A., P. Subramanian, K. P. Dere, and R. A. Howard 
(2000), Large- Angle Spectrometric Coronagraph Measure- 
ments of the Energetics of Coronal Mass Ejections, Astrophys. 
J., , 534 , 456-467, doi: 10. 1086/308747. 

Vourlidas, A., R. A. Howard, E. Esfandiari, S. Patsourakos, 
S. Yashiro, and G. Michalek (2010), Comprehensive Analysis 
of Coronal Mass Ejection Mass and Energy Properties Over a 
Full Solar Cycle, Astrophys. J., , 722, 1522, doi:10. 1088/0004- 
637X/722/2/1522. 



X - 12 



SAVANI ET AL.: MOMENTUM FLUX OF CMES 




Figure 1. Location of spacecraft during CME propaga- 
tion on the 27 March 2011. The shaded regions locate 
the field of view for the HI-1 cameras on both the ST- A 
and ST-B spacecraft. 
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Figure 2. A mass image from the COR2, HI-1 and HI- 
2 cameras from ST- A. The pixel intensity displays the 
locations of highest line of sight mass measurements. The 
cone overplotted on the figure displays the relevant part 
of the CME that propagates over the Earth. 





Figure 3. (Top) Panels show the graduated cylindrical 
shell (GCS) model overlaid onto COR2A and B mass 
images at 15.24, 25 March 2011. (Bottom) Shows the 
measurements of the CME mass in COR2 by using the 
GCS model as the CME boundary. 
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Figure 4. J-maps created for the position angle of the 
Earth (PA=97°). (Top) J- map displays the results from 
total brightness images with the location of Earth and 
the manually selected track of the CME's leading edge. 
(Bottom) Displays a mass J-map created from a sequence 
of mass images. 
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Figure 5. (Top)The time series of the WIND proton 
density and velocity for the period 25 March to 01 April 
2011. A normalised mass time series from the mass im- 
ages are overlaid onto the density axis (green). The mass 
time series is generated from the slit analysis at 3.8° elon- 
gation within the COR2-A camera. The time series is 
ballistically propagated to LI at a speed of 400km/s at 
the leading edge and 370km/s at the trailing edge, with a 
linearly decreasing speed profile. Qualitatively, the nor- 
mailsed mass time series shows a remarkable resemblance 
to the in situ density of the CME sheath. (Bottom) The 
mass J-map with the manually selected CME leading 
edge. 
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Figure 6. In situ data from ACE and WIND at LI point. 
The top 6 panels show the magnetic field in RTN carte- 
sian and spherical coordinates system; followed by proton 
velocity,density and Temperature, respectively. The ver- 
tical lines from left to right display the positions of the 
CME shock, and magnetic flux rope leading and trailing 
edge. 
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Figure 7. A comparison between the in situ data from 
LI (Bz component, velocity and density) to the averaged 
global response by the Dst index. The normalised mass 
time series from the COR2-A is displayed on the density 
axis as in figure 5. 
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Model at CCMC: BATSRUS 

Figure 8. 2D cut of the Earth's magneto-sphere and 
bow shock from the BAT-R-US simulation at CCMC. 
The colour scale of the density is the same in both panels. 
(Top)Displays the shape and structure of the magneto- 
sphere prior to the arrival of the CME shock. (Bottom) 
Shows the compressed magnetosphere at the time when 
the CME sheath is affecting the bow shock. 
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Figure 9. The time-derivative of the magnetic field mea- 
sured at the low-latitude Vassouras (Brazil) and Kakadu 
(Australia) ground stations. These stations represent the 
local noon and midnight of the CME shock front onto the 
bow shock respectively. The purple curves overlaid onto 
the data are the estimated values at these locations from 
the CCMC simulations. 




Figure 10. a) Shows the parameters from Thernisien 
et al. [2006] that are used to estimate the length of the 
CME along the line of sight, b) Is the initial volume 
element defined to be observed with the COR 2A camera, 
c) Shows the final dimensions of the volume element at 
LI, where the in situ measurements are taken. 



